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We analyse the extraction of the bound-state form factor from vacuum-to-hadron correlator, 
which is the basic object for the calculation of hadron form factors in the method of light-cone 
sum rules in QCD. We study this correlator in quantum mechanics, calculate it exactly, and derive 
the corresponding OPE. We then apply the standard procedures of QCD sum rules to isolate the 
ground-state form factor from this correlator. We demonstrate that fixing the effective continuum 
threshold, one of the key ingredients of the sum-rule calculation of bound-state parameters, poses a 
serious problem for sum rules based on vacuum-to-hadron correlators. 
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1. INTRODUCTION 

A calculation of hadron parameters from sum rules in QCD [l[ involves two steps: (i) one calculates the operator 
product expansion (OPE) series for a relevant correlator and obtains the sum rule which relates this OPE to the sum 
over hadronic states, and (ii) one attempts to get rid of the contribution of the excited states to the sum and isolate 
the ground-state contribution, thus relating parameters of the ground state to parameters of QCD. 

Constructing OPE in QCD is a well-defined procedure, whereas the extraction of the ground-state contribution 
requires several assumptions which do not come from the underlying field theory; it is hard to control the reliability 
of these assumptions and, as the consequence, the accuracy of the ground-state parameters extracted from sum rules. 

In order to probe the reliability and the accuracy of standard procedures adopted in the method of sum rules 
in QCD, quantum-mechanical model is an ideal tool: 1 in this model one can exactly calculate direct analogues of 
the field-theory correlators, one can generate analogues of the OPE for these correlators to any order, and, finally, 
one knows the exact bound-state parameters, such as masses, wave functions, and form factors. Therefore, one can 
compare the outcome of the standard sum-rule calculation with the exact values and in this way obtain an "unbiased" 
probe of the reliability of the sum-rule calculation. 

In the recent publications 0, Q i the procedures of extracting ground-state parameters from two-point [l[ and three- 
point 0] vacuum correlators were studied. In the present paper, a similar strategy is applied to vacuum-to-hadron 
correlator Q which is the basic object for the extraction of hadron form factors in the method of light-cone sum rules 
in QCD. 

We are going to demonstrate that fixing the effective continuum threshold for the case of vacuum-to-hadron corre- 
lator poses a challenging problem (in fact more challenging problem than that for the case of the vacuum-to- vacuum 
correlators). In particular, the effective continuum threshold is found to differ sizeably from nearly constant effec- 
tive continuum threshold for polarization operator and to depend strongly both on the momentum transfer and on 
the Borel parameter. This leads to difficulties in finding a reasonable approximation to this quantity and, as the 
consequence, to large systematic uncertainties in sum-rule results for hadron form factors. 

2. THE MODEL 

We consider a non-relativistic model with HO potential 

H = H + V(r), H = f/2m, V(r) = mwV/2, r = \r\. (2.1) 
The full Green function G(E) = (H — E)^ 1 and the free Green function Go(E) = (Ho — E)^ 1 satisfy the equation 

G- 1 {E)-Go\E)=V 1 (2.2) 



1 For a discussion of many aspects of sum rules in quantum mechanics we refer to 6, 7, 8, 9]. 



2 



which may be solved by constructing the expansion in powers of the interaction V: 

G{E) = G (E) - G (E)VG (E) + ■■■ . (2.3) 

In HO model all characteristics of the bound states are easily calculable, e.g. for the ground state one has 

3 /rnaj\ 3 / 4 ( \ \ 

Aground = x^) *ground(^) = \—) ex P ( ~ -jiiiuur 2 J , F ground (q) = cxp(-q 2 /4mw) . (2.4) 

The elastic form factor of the ground state F gmun( ±(q) is defined according to 

F gIoun d{q) = (Aground | | Aground) = / *^ round (fc)V> gr ound(fc - q) = J dr\4! gIound (r)\ 2 e lqr , q=\q\ (2.5) 

with the current operator given by the kernel 

{f\J{q)\f) = exp(iqr)6{r - f). (2.6) 

3. THE POLARIZATION OPERATOR 

The polarization operator 

n(T) - (f f = 0| exp(-HT)\f = 0} (3.7) 

is used in the method of sum rules for the extraction of the wave function at the origin (the decay constant) of the 
ground state A detailed discussion of the corresponding procedure for HO model was presented in For HO 
potential, the analytic expression for n(T) is known [6j: 

n(T) = m 3/ ' (3-8) 



V 7T / 



[2 sinh(u;T)] 3/2 



The corresponding average energy is defined as follows 



Ejx(T) = -a T logn(T) = ^coth(wT). (3.9) 
The OPE is the expansion of the polarization operator at small Euclidean time T (or, equivalently, in powers of ui): 

3/2 
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(3.10) 



THE VERTEX FUNCTION 



The basic quantity for the extraction of the form factor in the method of dispersive sum rules is the correlator of 
three currents [J] . The analogue of this quantity in quantum mechanics has the form 

rfo.Ti,?) = (f f = 0\exp(-Hr 2 )J(q) exp{~H Tl )\n = 0), (4.11) 

with J(q) defined in (|2.6[) . For equal times t\ = t 2 = T/2, we obtained an explicit expression for HO model 0] 

T(T,q) = « 3/2 1 exp L^tanh(cT/2)) . (4.12) 

v k ' [2sinh(wT)] J/i V 4mw / 

The corresponding average energy has a simple form 

E r (T,q) ee -d T \ogT(T,q) = Lcoth(a;r) + j- - 1 (4.13) 

2 4m (1 + cosh(cjJ )) 

The OPE is the expansion of T(T, q) in powers of o>; for its explicit form in HO model we refer to [3|]. An unpleasant 
feature of Tope is the polynomial growth of power corrections with q 2 which limits the convergence of the truncated 
Tope to not very large q 2 . Respectively, this correlator may be used for extracting form factors at not very large q 2 . 
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5. THE VACUUM-TO-HADRON CORRELATOR 



To apply the method of sum rules to hadron form factors in a wider range of q 2 , a different correlator is used - the 
vacuum-to-hadron amplitude of the T-product of two currents. The analogue of this quantity in quantum mechanics 
has the form 

A(E,q) = (r = 0\G(E)J(q)\V gIOUIld ), (5.14) 

or, after Borelization, 

A(T,q) = (r = 0\G{T) J(g)|* ground ). (5.15) 

An obvious disadvantage of this correlator compared to the three-point function T is the necessity to know the ground- 
state wave function for its calculation. As we shall see, as a bonus, this correlator has an enhanced contribution of 
the ground state which makes it potentially more attractive for the extraction of the ground-state form factor. 

In HO model, the ground state wave function is known explicitly. Making use of this wave function (|2.4p and the 
results for the Green function in configuration space G(r,T) = (f = 0\G(T)\r) from 0], we obtain the exact analytic 
expression for A in HO model: 

«^(7)'"-(-H-(-ii"-'i' < 5 - ie > 

The function A(T, q) depends on two dimensionless variables q 2 /Amui and ojT. The average energy for this correlator 
has a simple form 

3 a 2 

E(T,q) = -St log A(T,q) = -u + ±- exp(-2wT). (5.17) 

2 Zm 

It should be noticed that the properties of the correlator A(T, q) differ considerably from the properties of n(T) and 
F(T, q): for instance, A(T,q) and the corresponding E(T,q) are finite at T = 0, whereas the correlators n and T 
diverge at small T as T~ 3 / 2 . The energies, related to n and T, Eji and Er, diverge as 1/T. 



A. Ground-state contribution to A(T, q) 



The ground-state contribution to the correlator reads 

Aground(T, q) = Aground (f = 0) CXp (--EgroundT) F groun d(q) (5.18) 

leading to 



Aground (T,q) ( q 



-,2 



exp I — ^—e- 2 ^ T . (5.19) 
A(T, q) y \ Amu; J V ' 

Notice the following features: 

(i) At large T, the ground state provides the dominant contribution to A, similar to any other correlator: 

A(T, q) -f * ground (r - 0)e- £ — dT F ground (q) + . . . (5.20) 

(ii) At small q, the ground state dominates the correlator for all T. This is a specific feature of A which arises due 
to the choice of the initial state. Thus, compared with H(T) and T(T,q), the correlator A "maximizes" the ground 
state contribution. 



B. OPE for A(T, q) 

In order to obtain the OPE for A along the lines of QCD, one should expand the Green function G in Eq. (|5. 15[) in 
powers of the interaction (Fig. [l}. This generates the expansion of A(T,q) in powers of (uiT) 2 . To obtain this series 
from the exact expression requires however some care: one should take into account that the parameter u> enters both 
the wave function and the Green function G. Therefore, in practice it is convenient to proceed as follows: treat the 
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A (E,q) A 2 (E,q) A 4 (E,q) 

Fig. 1: Expansion of the Green function G(E) in the correlator A(E,q) in powers of the interaction. Full blobs denote the 
wave function of the ground state, dashed lines correspond to the potential V. 

parameter in G as u, but the same parameter in the wave function as luq] calculate the exact correlator A(T, q\u>, u>o) 
from Eq. (15. 15|) : expand this expression in powers of u> and obtain A(T, q\u>, ujq) = 'Y^ = ${uT) 2n A2 n {T , q\jJo)\ finally, 
set ojo —> w in the functions A2 n (T, q\u>o)- This procedure yields a rather complicated OPE series: 

oo 

A(T,q) = J2("T) 2n A 2n (T,q\oj), (5.21) 

n=0 

where the functions A n (T, q\u>) have a nontrivial T, q, and w-dependence: 
/mw\3/4 1 / q 2 uT 



7T / (1 + wT) 3 / 2 V 2muj 1 + 



12 (-) JTT^W 2 exp (-^TT^ j t 3(1 + wT)(3 + wT) - 2q T l m \ ■ (5-22) 



For large q , one obtains 



'■-'■'''""'-'■(-^il^l 



Obviously, for large g 2 all functions A„ remain finite and tend to zero, so that the truncated OPE may be calculated 
for any value of q 2 . However, A n+2 /A n ~ q 2 and thus, for a purely confined potential, the truncated OPE does not 
reproduce the exact A(T,q) in the limit q 2 — > oo. 2 Fig. [2] illustrates the accuracy of the OPE series, truncated to 4 
terms, Aope.4 = ^2, n =a{uT) 2n A2 n . In order to have a good description of the exact function A(T, q) for q < 8ui with 
^4ope,4 one should stay in the region uiT < 1 ^ 1.2. 




Fig. 2: The accuracy of the OPE series truncated to 4 terms Aope,4 = ELo( w ^) 2 " A 2n f° r A(T, q) vs q for two different values 
of luT: (a) cuT = 1.0, (b) ljT = 1.5. Dashed (red) line - A , dotted (blue) - A 2 {ujT) 2 , long-dashed (green) - A^loT) 4 , solid 
(violet) - A 6 (u/T) 6 , Thick solid (black) line - 4qpe,4. 



2 Compare with the following example: imagine we expand / = e~ 2x = e~ x e~ x = (1 — x + a; 2 /2 + . . .)e~ x = JZ^Lo fn( x )- For x — + oo, 
fn(x) —* 0, so any truncated series takes a finite value for x — » oo. However, /n+l (%)/ fn (x) — > oo, and thus a truncated series does not 
reproduce the exact f(x) for large x. 
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Zeff (gT = 0-5, q) z eff (o)T= 1, q) 




Fig. 3: The exact effective continuum threshold z e s(T, q) [as obtained by solving (15.28[) for the exact bound-state parameters 
in the l.h.s.] vs q for uiT = 0.5 (a) and uiT = 1 (b). A horizontal line at 2.5 is given as a benchmark: this constant was shown 
to provide a good approximation to the exact continuum threshold in sum rule for H{T) [|. 



C. Sum rule and exact effective continuum threshold 



The sum rule represents the equality of the correlator calculated in quark and in hadron basis: 

A (T, q) + A powcl (T, q) = A gIound (q) + A cxcitcd (T, q). (5-24) 

According to the standard procedures of the method of sum rules, the contribution of the ground state is assumed to 
be dual to the low-energy region, i.e. the region of small momentum transfers. To implement this idea technically, one 
should represent the correlator as the dispersion representation in z = fc 2 , and then cut the integral at the effective 
continuum threshold z c g. The easiest way to derive the spectral representation for A{T 1 q) is to use the full Green 
function in momentum space 

G(i?,T) s (f =0| G (T)|S> = {2w ) sl2[co3h( l T)]sl ^ '-^ (5-25) 

The correlator A(T, q) takes the form 

A(T, q)= fdk G(k 2 , T)* ground (fc - q) = / dze~ zT a{z, q) (5.26) 



The explicit expression for a(z, q) may be easily obtained using Aground given in 
The cut correlator is denned according to 



A cnt (T,q,z cfi )= e- ZI a(z,q)dz, (5.27) 
Jo 

and it should reproduce the ground-state contribution: 

VWnd(r = 0)e- £ — dT F ground (q) = A cut (T, q, z cS (T, q)). (5.28) 

This expression is precise if we use an exact T- and g-dependent effective continuum threshold (which cannot be 
calculated knowing only OPE, but can of course be reconstructed in HO model since we know the exact form factor, 
ground-state energy, and ground-state mass). Moreover, Eq. (|5.28|) may be understood as the definition of the exact 
effective continuum threshold, if one makes use of the exact hadron parameters on the l.h.s. Obviously, Eq. (|5.28p 
alone is not predictive and the form factor (as well as any other parameter) of the ground state may be obtained 
in the method of sum rules only if one imposes an independent criterion to fix the effective continuum threshold. 
Recall that this criterion does not come directly from the underlying theory. One may expect that in the case of 
the correlator A(T, q), the exact z e g strongly depends on the momentum transfer: e.g., as noticed above, the ground 
state dominates the correlator at small q, therefore z c g(T, q) — > oo as q — > 0. Fig. §3§ shows z c g(T, q) as obtained by 
solving (|5.28p for the exact parameters of the ground state from (|2 .4[) . Indeed, the exact effective continuum threshold 
strongly depends both on T and q. Moreover, it is rather far from the value z c = 2.5 u which was found to provide a 
good approximation for the exact continuum threshold for two-point sum rule Q ■ 

Recall that the effective continuum threshold is one of the two basic ingredients (along with the OPE for the 
correlator) which determine the value of the ground state parameter. Therefore, the observed strong dependence of 
the exact effective continuum threshold on q and T is very unpleasant: it implies difficulties in finding a reasonable 
approximation to this quantity and, as the consequence, large uncertainties in the sum-rule results for the hadron 
form factor. 
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D. Extraction of the ground-state form factor 

Let us now apply the standard procedures and extract the form factor at two values of q: at the intermediate 
q~2u (q 2 /4muj = 1) and at relatively large q — 4w (q 2 /4mw = 4). 

In realistic cases, the exact correlator is not known and one may calculate only several terms of its OPE series. 
Assume we have 4 terms at our disposal: 

3 

^ope,4 = ^(u/T) 2 "A 2n (T, q). (5.29) 

First, let us determine the "window" - the region of the Borel parameter where, according to standard criteria, we 
may work to extract the ground-state form factor: To be sure that our truncate OPE provides a good description of 
the exact expression - say, with 1% accuracy - we must stay at not very large values of T. This requirement gives 
the upper boundary of the Borel "window", ljT < 1.2 The lower boundary of the "window" is determined by the 
requirement that the gound state gives a sizeable contribution - say, more than 50% - to the correlator. This yields 
0.25 < loT for q = 2lo and 0.8 < ujT for q = Auj. 

So the "window" is 0.25 < uoT < 1.2 for q = 2u and 0.8 < uT < 1.2 for q = 4w. 3 

Second, we need to apply the cut to the correlator. Let us define the cut correlator in a bit different way than in 
the previous Section, to make it closer to the standard sum-rule calculation in QCD: namely, we apply the cut to Aq 
and leave higher terms intact. Then the cut correlator has the form 

A mt (T,q,z c ) = A c ut (T,q,z c )+A 2 (T,q)+A 4 (T,q)+A e (T,q), Ay nt (T, q, z c ) = f " dz exp(-zT)a (z, q). 

Jo 

(5.30) 

Here ao is the spectral density of Aq: 

/r°° i 
dkG (k 2 ,T)V smund (k-q) = J dze- zT a a (z,q), G (k 2 , T) = ^-^ e ~ Tk2/2m . (5.31) 

The form factor is related to the cut correlator according to relation l|5.28[) . 

Next, the key step - to impose a criterion for fixing z c . In the literature, two prescriptions for fixing z c have been 
used: 

(A) One assumes z c to be a T- and g-independent constant and its value is chosen to be the same as in the two-point 
sum rule [lol ]. In our case this implies setting z c = 2.5 u. 

(B) One assumes z c to be T-independent, but tunes its value for each value of q, separately. To this end one calculates 
the average energy of the cut correlator 

E(T,q,z c ) = -d T log A(T,q,z c ), (5.32) 

which depends on T and q because of the approximation z c g(T,q) — > z c (see discussion in Q). Then, one determines 
z c for each q such that the function E(T,q, z c ) has a horizontal tangent E = -Eground [LH. 

We also present the exact effective continuum threshold, which is obtained by solving Eq. (]5 . 28[) for the known 
ground-state parameters on its l.h.s., and the correlator with the cut applied according to (|5.30[) on its r.h.s. 

The corresponding plots are given in Figs. 2] and [5] According to the usual procedures, the theoretical uncertainty 
of the sum-rule prediction is obtained from the range covered by the form factor when varying the Borel parameter 
within the "window" . 

For q = 2uj (Fig. the exact effective threshold turns out to depend sizeably on T in the "window" . Prescription 
B leads to a stable form factor -FsR/^ground = 0.82 ± 0.02 which however underestimates the actual value by almost 
20%; prescription A leads in this case to a very good estimate Fsr/ '-^ground = 1-02 ± 0.03. 

For q — 4uj (Fig. [5]) the situation is different: The exact continuum threshold is almost flat in the "window", however 
its value is different from the values determined by both prescriptions A and B. As the result, the prescription A 
strongly underestimates the form factor, leading to -FsR/fground = 0.91±0.01, whereas the prescription B overestimates 
it, yielding iW-Fground = 1.08 ± 0.02. 



3 Notice that the "window" disappears for large values of q: with a fixed number of terms in the truncated OPE, the upper boundary of 
the "window" decreases with q, whereas the lower boundary increases, as can be seen from (15.191 1. 
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Fig. 4: Intermediate momentum transfer q — 2lj (q 2 /Amu) = 1). Dashed (blue) line - constant effective continuum threshold 
z c — 2.5 lu (prescription A); dotted (red) - T-independent effective continuum threshold tuned according to prescription B; solid 
(black) - T-dependent exact effective continuum threshold z c s(T,q). 



z ef f (T, q) 



3.4 
3.2 
3 
2.8 
2.6 
2.4 
2.2 



0.2 0.4 0.6 0.8 



a/I 



E(T, q, Zc) 



2.5 
2.25 
2 

1.75 
1.5 
1.25 
1 

0.75 



/\ 








1 \ 








\\ 








\ „ - - 




— ' 





0.2 0.4 0.6 O.i 



ml 



t-'sn IT. q. z L .) 
F„ (q ) 



1.1 
1 

0.9 
0.8 
0.7 







\ 





0.2 0.4 0.6 0.8 1 



ml 



Fig. 5: Relatively large momentum transfer q = Ad (q 2 /Amu — 4). Dashed (blue) line - constant effective continuum threshold 
z c — 2.5u (prescription A); dotted (red) - T-independent effective continuum threshold tuned according to prescription B; solid 
(black) - T-dependent exact effective continuum threshold z c ff(T,q). 



The lesson to be learnt from these two examples is clear - neither the recepy A or B can guarantee the extraction 
of the true form factor with a controlled accuracy: in "good" cases one of these prescriptions may work, in "bad" 
cases they fail. Unfortunately, so far we have not found any criterion which could discriminate between the "good" 
and the "bad" cases. 



6. CONCLUSIONS 



We studied the vacuum-to-hadron correlator A(T, q) in harmonic-oscillator model and analysed the procedure of 
extracting the ground-state form factor from this correlator using the procedures of the method of sum rules. The 
great advantage of this model is the possibility to obtain exact expression for the correlator and the knowledge of 
exact parameters of the ground state. Therefore this model provides a brilliant opportunity to study the reliability of 
the procedures used in the method of sum rules. 

Our main conclusions are summarized below: 

1. Compared to vacuum-to- vacuum correlator, the correlator A(T, q) requires the knowledge of the ground-state wave 
function for its calculation. As a bonus, it contains an enhanced contribution of the ground state, which makes it 
potentially a favourable candidate for the calculation of the ground state form factor with sum rules. 

2. However, the extraction of the ground-state form factor from A(T,q) faces a challenging problem of fixing the 
effective continuum threshold: as demonstrated in our exactly solvable model, the exact effective continuum threshold 
- a solution of the sum rule for the known exact ground-state parameters - strongly depends both on the Borel 
parameter T and the momentum transfer q. This implies serious difficulties in finding a reasonable approximation to 
this quantity and, as the consequence, large systematic uncertainties in sum-rule results for hadron form factors. 

3. We tested prescriptions used in the literature for fixing the effective continuum threshold and found that none of 
these prescriptions guarantees the reliable extraction of the form factor: for some values of the momentum transfer 
they may provide a good result, for other values of the momentum transfer they may fail. Unfortunately, so far we 
have not found any criterion which could discriminate between the "good" and the "bad" situations. 

4. We have seen that the usual criteria adopted in the method of sum rules for obtaining error estimates of the ground- 
state form factors (Borel stability, variation of the form factor in the "window" ) do not guarantee a good extraction 
of the form factor and do not provide realistic error estimates - the actual errors turn out to be considerably larger. 
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In this respect hadron form factors obtained from the vacuum-to-hadron correlators share the same difficulties as 
hadron parameters obtained from other versions of sum rules (see conclusions in @, 9)- 

A systematic study of the actual accuracy of hadron form factors obtained from vacuum-to-hadron correlators (in 
particular, with light-cone sum rules) should be undertaken before these form factors may be used in problems where 
rigorous error estimates are required. 
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